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We review some recent results for a system of simple random walks on graphs, known as frog model. Also, we discuss 
several modifications of this model, and present a few open problems. A simple version of the frog model can be 
described as follows: There are active and sleeping particles living on some graph. Each active particle performs a 
simple random walk with discrete time and at each moment it may disappear with probability 1 — p. When an active 
particle hits a sleeping particle, the latter becomes active. 

Keywords: simple random walk, critical probability, shape theorem, recurrence 


1 Introduction 

The purpose of this note is to describe recent results and to propose some problems for future research 
for the system of interacting random walk called frog model. The basic version of this model can be 
informally described in the following way. Initially there is a random number of particles at each site of a 
graph Cj. A site of (j is singled out and called its root. All particles are sleeping at time zero, except for 
those that might be placed at the root, which are active. At each instant of time, each active particle may 
die with probability (1 — p). Once an active particle survives, it jumps on one of its nearest neighbors, 
chosen with uniform probability, performing a discrete time simple symmetric random walk (SRW) on 
(j. Up to the time it dies, it activates all sleeping particles it hits along its way. From the moment they 
are activated on, every such particle starts to walk, performing exactly the same dynamics, independent 
of everything else. 

In this paper we consider several modifications of this model always keeping the following two common 
features: 

• every particle performs a simple random walk, independently of everything else (although the mo¬ 
ment when it starts may be random); 

• there is no influx of particles to the system (but the particles may die), which distinguishes the 
models we consider here from e.g. branching random walks. 

^The author is thankful to CNPq (302981/02-0) for partial support. 
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The model described above is a modification of a model for information spreading first proposed by 
K. Ravishankar. The idea is that every active particle has some information and it shares that information 
with a sleeping particle at the time the former jumps on the latter. Particles that have the information move 
freely helping in the process of spreading the information. Combustion chemical reactions provide another 
motivation for studying this process, cf. [ 8 ]. The model that we deal with in this paper is (with p = 1) a 
discrete-time version of that proposed by R. Durrett (1996, private communication to F. Machado), who 
also suggested the term “frog model”. Our study of the frog model with p < 1 originated from a question 
asked by I. Benjamini to F. Machado in 2000. 

To the best of our knowledge, the first published result on this model (with p = 1, Cj = I' 1 and one- 
particle-per-site initial configuration) is due to Teles and Wormald [9] where it was referred to as the “egg 
model”. They proved that, starting from that initial configuration, the origin will be visited infinitely often 
a.s. See Section 4 for more results about the probability of visiting the origin infinitely often. 

Now, we give a formal definition of the model described above. 

Let q be a random variable taking values in N = {0,1,2,...} such that P[tj > 1] > 0. Let {q(x);x G 
£},{(S*(/))„ eN ;/ G {1,2,3,...},x G Cj} and {(E*(/));/ G {1,2,3,...},* G £ 7 } be independent sets of 
i.i.d. random variables defined as follows. For each x G Cj, q(x) is an independent copy of q, and gives 
the initial number of particles at site x. If q(x) > 1, then for each 1 < i < q(x), (S*(i))„ £ N is a discrete 
time SRW on Cj starting from x (it describes the trajectory of /-th particle departing from x), and E*(/), 
which stands for the lifetime of /-th particle departing from x, is a random variable whose law is given 
by P[Ep(/) = k] = (1 — p)p k ~ l , k= 1,2,..., where p G [0,1] is a fixed parameter. Thus, the /-th particle 
at site x follows the SRW (S£(/)) ne n and dies (disappears) S x (i) units of time after being activated. For 
x ^ v let 

t(x,y) = min min{« < S*(i) : S x n (i) = y} 

1<Kt|(.y) 

(note that, if the SRW is transient on Cj, then t(x.y) = °° with positive probability). The moment when all 
the particles in x are awakened is defined as 

r( 0 ,x) = inf{f(x 0 ;*i) H-hf(x m _t,x m )}, 

where the infimum is over all finite sequences 0 = xo,xi,... ,x m _i ,x m = x. Observe that T (0,x) = °° means 
that the site x is never visited by active particles. 

It is important to note that at the moment the particle disappears, it is not able to activate other particles 
(as first we decide whether the particle survives, and only after that the particle that survived is allowed to 
jump). Notice that there is no interaction between active particles, which means that each active particle 
moves independently of everything else. We denote by FM(^,/j,q) the frog model on the graph Cj with 
survival parameter p and initial configuration given by independent copies of q in each site of Cj. 

Next, we review known results for frogs, and propose some open problems. 

2 Shape theorem 

In this section we consider only the case p = 1, i.e., the particles never die. 

Let be the set of sites which were visited by active particles up to time n, namely, = {y G Z d : 
T(0,y) < n}. Define also \ n = {>’ + (- l/2,\/2] d :>• G q„} C W 1 . 
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Theorem 2.1 ([3], Theorem 1.1). Consider the model FM(Z‘ I , 1,1]). For any dimension d > 1 there is 
a nonempty convex set FL = FL(d, T|) C W 1 such that for almost all initial configurations, conditioned on 
|r|( 0 ) > 1 }, we have for any 0 < £ < 1 


P 

(1 -e)FLc — c (l + e)FL 

n 


for all n large enough a.s. 

For the particular case r| = 1 this result was obtained in [1] (Theorem 1.1). Naturally, the main tool 
for proving the above shape result is the Subadditive Ergodic Theorem. When applying that theorem, the 
most difficult part (especially in the dimension d > 2) was to prove that the mean time until some active 
particle comes to a site xo is finite for all xo- 

Note that, although Theorem 2.1 establishes the existence of the asymptotic shape FL, it is difficult to 
identify this shape exactly. Of course, FL is symmetric, the origin belongs to the interior of FL, and FI C ‘D, 

where T) = {x = (x^,... ,x^) g : |x^| -\ -h |x^| < 1). Also, note that if the initial configuration 

is augmented (i.e. some new particles are added), then the asymptotic shape (when it exists) augments as 
well. We are going to show that if the initial configuration is heavy enough, then the limiting shape FI 
may contain some pieces of the boundary of ‘D (a “flat edge” result), or even coincide with ‘D (a “full 
diamond” result). 

To formulate the “flat edge” result, we need some additional notation. For d >2 and 1 < i < j < d let 
A ij = {x = (x«,... jX^) e W l : xW = 0 for k ± i,j}, and for 0 < p < 1/2 let 

©f- = {x = (x^,... ,xW) e A ij : |x«| + |xW| = l,min{|x^|, |x^|) > p}. 


Define ©P to be the convex hull of (©f 7 )i</<;<d- We are interested now in the asymptotic shape in the 
frog model when the initial configuration is such that any site xgS ' 1 contains exactly m particles. 

Theorem 2.2 ([1], Theorem 1.2). Consider the model FM(Z el , 1 ,m). For each d >2 there exists mo = 
mo(d) such that for all m > mo there exists 0 < p < 1/2 such that ©P C FL(d,m). 

It is not difficult to see that, by monotonicity, the above result is verified for FM(Z‘ I , l,r|) with r| such 
that r| > mo a.s. 

Next, we formulate the “full diamond” result: 


Theorem 2.3 ([3], Theorem 1.2). Consider the model FM(Z d , 1,1]). Suppose that for some positive P < d 
and for all n large enough we have 


P[t](x) > n] > 


1 

(logn)P' 


Then, Theorem 2.1 is verified with FL = T). 

Remark. In [ 8 ], a continuous-time analogue of FM(Z el , 1,1) was considered. They prove the shape 
theorem (using the method different from that of [1, 3]), and the fact that the local particle distribution 
tends to the product of Poissons. Also, it is proven there that when the dimension is large enough, the 
asymptotic shape is not the Euclidean ball (the speed in the axial direction is greater than the speed in the 
diagonal direction). 
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3 Phase transition 

For the model FM(( 7 ,/?,r|), let us consider the following definition. 

Definition 3.1. A particular realization of the frog model survives if for every instant of time there is at 
least one active particle. Otherwise, we say that it dies out. 

Now we observe that P[FM(( 7 ,/?,r|) survives] is nondecreasing in p and define 

Pc{Q, r|) := inf{p : P[FM(£,/?,r|) survives] > 0}. 

As usual, we say that FM(( 7 ,p,r|) exhibits phase transition if 0 < p c (ff r|) < 1- 

Now we present two lower bounds on p c (tf,r |) which are obtained by a direct comparison with a 
Gabon-Watson branching process. 

Proposition 3.1 ([2], Propositions 1.1 and 1.2). If Ft] < °o then for arbitrary graph (f it holds that 
Pc{Q, T|) > (Er| + 1) _ . If one makes an additional assumption that (f is a graph of maximum degree k, 
then 

~ l + (Ar—l)(Er| + l)' 

Apart from comparison with a Gabon-Watson branching processes, the main idea for proving most of 
the results of Section 3 is to think of the frog model as of ( dependent) oriented percolation model. Namely, 
for each x, we consider its “virtual range”, i.e., the set of sites visited by the particles originating from x, 
provided that they are eventually activated. To construct the oriented percolation process, we then draw 
oriented edges from x to all sites belonging to the virtual range of x. 

3.1 Extinction and survival 

Next, we present some more general results concerning the extinction and the survival of the process, 
mainly for the cases of 7L d and regular trees. 

We begin by showing that, under mild conditions on the initial number of particles, the process dies out 
a.s. in Z for every p < 1. From now on, a V b stands for ma x{a,b}. 

Theorem 3.1 ([2], Theorem 1.1). 7/Elog(r| V 1) < then p c ( Z,r|) = 1. 

Next, we find sufficient conditions to guarantee that the process becomes extinct for p small enough in 
YL d , d > 2, and in T f /, d > 3 (by we mean a tree such that the degree of each vertex is d). Namely, for 
the case of tree, the result is given by 

Theorem 3.2 ([2], Theorem 1.2). Suppose that there exists 8 > 0 such that Er| 8 < °°. Then /?, (T,/. r|) > 0, 
i.e., the process on dies out a.s. for p > 0 small enough. 

For the case of ^-dimensional lattice, the condition for the positivity of p c ( Z d ,r\) is described by 
Theorem 3.3 ([2], Theorem 1.3). Suppose that E(log(r| V l)) f/ < °°. Then p c (Z d ,r\) > 0. 

Now, let us state the results related to the survival of the process. The next theorem shows that for any 
nontrivial T| the frog model survives on Z d , d > 2, and on T,/, d > 3, when the parameter p is close enough 
to 1. 

Theorem 3.4 ([2], Theorems 1.4 and 1.5). If P[r| > 1] > 0, then p c {Z d ,r\) < 1 for all d >2 and 
p c (T,/,r|) < l for all d> 3. 
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The method of proving Theorem 3.4 is based, of course, on the fact that the the usual percolation models 
on Z d , d > 2, and on T,/, d > 3, percolate when the parameter is close enough to 1. 

The next result is the counterpart of Theorem 3.2. Note that Theorems 3.2 and 3.5 give the complete 
classification in r| of the frog model on T)/ from the point of view of positivity of p c (T<i,r|). 

Theorem 3.5 ([2], Theorem 1.6). //Er| 8 = °° for every 8 > 0, then /?, (T,/. r|) = 0. 

Besides, we are able to show that, for every fixed d, if r) has a sufficiently heavy tail, then FM(Z rf ./;. r|) 
survives with positive probability for all values of p € (0,1) (which would be the counterpart of Theo¬ 
rem 3.3). However, we do not state this result now, since we will give a stronger result (cf. Theorem 4.3) 
later. 

3.2 Asymptotic values for the critical parameter and monotonicity 

The following theorem gives asymptotic values for critical parameters (compare with Proposition 3.1) for 
the case of Z d and regular trees. 

Theorem 3.6 ([2], Theorems 1.7 and 1.8). We have, for the case Er| < °° 

lim p c ( T rf ,r|) = lim p c (Z d , r|) = 1 ■ 

cf-joo d-toc i + H,T| 

Remark. Observe that by truncating T| and using a simple coupling argument one gets that if Er| = °°, then 

lim /^.(T^ri) = lim p c (Z d ,r |) = 0. 

d-f oo d—} °o 


Note that Theorem 3.6 together with Proposition 3.1 suggest that there is some monotonicity of the 
critical probability in dimension. Then, a natural question to ask is the following: 

Open Problem 3.1. Is it true that p c (Z d ;r |) > p c (Z d+l ,r|) for all d (or at least for d large enough; also, 
can one substitute “>” by “> ”)? What about the case ofTj? 

In fact, there is a more general question: for Q\ C ( 72 , is it true that p c (Qi,r[) > p c (<f 2 ,t\), as hap¬ 
pens for usual percolation models? In a recent paper [5] that last question was answered negatively. To 
formulate their result, let us call (f*(n) the graph whose set of vertices is U A . €Z 2 N n (x) where N„(x) = 
{(x,0), (x, 1),.,., (x,n)}. Moreover, its set of edges is such that for any x £ Z 2 the vertex (x. i) is neighbor 
of (x, j) for 0 < i < j <n. Besides, (x,0) is neighbor of (y,0) if x.y are neighbors in Z 2 . In words, to 
construct the graph <f*(n), one have to attach complete graphs with n vertices to each site of Z 2 . 

Theorem 3.7 (see [5]). There exists no such that for all n > no we have p c (Z 2 , 1) < p c (Cj*(n), 1). 

4 Recurrence and transience 

Again, consider the model FM (Q,p. r\). One may also be interested in the study of other types of critical 
behaviour with respect to the parameter p. Consider the following 

Definition 4.1. The model FM( (j ,p. x\) is called recurrent if 

P [the root (the origin) is hit infinitely often in FM(£ 7 ,/ 2 ,T|)] > 0. 

Otherwise, the model is called transient. 
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Note that, even in the case when a single SRW on Q is transient, it is still reasonable to expect that 
the frog model with p = 1 is recurrent. For example, for the model FM(Z d , 1,1) the recurrence was 
established in [9]. However, establishing the recurrence property for that model is nontrivial; for example, 
consider the following 

Open Problem 4.1. Is the model FM(T,/. 1,1) recurrent? 

Now, denote 


p„(f, r|) = inf{p : P[the root is hit infinitely often in FM((7,p,r|)] > 0} 


(here, by definition, inf0 = 1); clearly, p u (Cj .r\) > p c { (j.r\) for every (f and rp Now, we are interested in 
studying the existence of phase transition with respect to p u . 

First, we discuss some situations when the model is transient for every p, except possibly the case 
P = I- 

Theorem 4.1 ([2], Theorems 1.9 and 1.10). Suppose that Er| e < °° for every 0 < £ < 1. Then /;„(T,/, r|) = 
1. Now, suppose that E[log(r| V l)] f/ < °° Then p u (Z d ,r\) = 1. 

The next two theorems give sufficient conditions to have p u < 1 on trees and on Z d . 

Theorem 4.2 ([2], Theorem 1.11). Suppose that there exists p < such that 

p h > »] > 4 


for all n large enough. Then T|) < 1. 

Theorem 4.3 ([2], Theorem 1.12). Suppose that there exist p < d such that 


P[r| > n\ > 


1 

(logn)P 


for all n large enough. Then p u (Z d , r|) = 0. 

(Compare this result with Theorem 2.3.) 

Theorems 4.2 and 4.3 give sufficient conditions on the tail of the distribution of r| for the process to 
be recurrent when p < 1. On the other hand. Theorem 4.1 shows that for the one-particle-per-site initial 
configuration the process is not recurrent even if the parameter p < 1 is very close to 1. The model 
with one-particle-per-site initial configuration being the most natural example one has to hand, a natural 
question is raised: What can be done (i.e., how can one modify the model) to make the model recurrent 
without augmenting the initial configuration? Notice that, by definition, in our model the lifetime of 
active particles is geometrically distributed. In order to find answers to that question, we change this and 
consider the situation when the lifetime has another distribution, possibly a more heavy-tailed one. 

Let E be a nonnegative integer-valued random variable. Let us consider the frog model on (f with one- 
particle-per-site initial configuration, and the lifetimes of the particles after activation are i.i.d. random 
variables (E x ,x £ If) having the same law as E. This model is denoted by FM(( 7 ,E). 

Theorem 4.4 ([2], Theorem 1.13). Suppose that one of the following alternatives holds: 


f = Z and E\/E < 
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• $ = z2andE T^kr) <°°’ 

• = Td or Z cl , d > 3 and E E < 

Then FM(^,E) is transient. 

Theorem 4.5 ([2], Theorem 1.14). For each dimension d there exists (3^/ > 0 such that if for all n large 
enough one of the following alternatives holds 

• d = 1 and P[E > n 2 ] > pi« _1 log log n, 

• d = 2 and P[S > rf \ > P 2 « _2 (log«) 2 , 

• d > 3 and P [Z>n 2 ]>& d n 2 log n, 
then FM(Z d , E) is recurrent. 

Theorems 4.1 and 4.4 are proved by making all the particles initially active, and getting an upper bound 
on the mean number of sites which send a particle to the origin. Theorems 4.2, 4.3, and 4.5 also are 

proved using a common approach: We think of (f as a disjoint union of sets 3k-k = 1,2,_of increasing 

sizes, such that with large probability (increasing with k), the set 3k contains a lot of particles in the initial 
configuration. Besides, given that 3k contains many particles, also with large probability (increasing 
with k as well), the virtual paths of those particles will cover the whole set 3k+\ together with the origin, 
thus activating all particles placed originally in 3/.+ I • With a particular choice of that sequence of sets, 
the intersection of the events mentioned above occurs with strictly positive probability, which implies, 
consequently, that the process is recurrent (as in this case for each k there is a particle from 3k which visits 
the origin, and so the total number of particles visiting the origin is infinite). 

As noted earlier in this paper, the question of recurrence of frog models is nontrivial even for the case of 
infinite lifetime (when, of course, the random walk of individual particle is transient). In view of that, the 
following problem is of interest. Given a collection of numbers (0 < q(x) < 1 ,x £ Z d \ {0}), d > 3, put a 
sleeping particle into site x with probability q(x). The configuration of sleeping particles is then fixed, and 
the process is started. The main question for such a model is: How can one distinguish transience from 
recurrence by looking at the initial density q(x)2 The next theorem shows that the critical (i.e. separating 
transience from recurrence) rate of decay of the function q(x) is ||x|| -2 : 

Theorem 4.6 ([7], Theorem 1.1). There exists a cr = OC cr (d), 0 < 0C cr < °°, such that 

(i) if a < a C r and q(x) < a||x|| -2 /or all x large enough, then the process is transient; 

(ii) if a > (X C r and q(x) > a||x|| -2 /or all x large enough, then the process is recurrent. 

The first part of Theorem 4.6, is proved by domination by branching random walk; the proof of the 
recurrence part is, however, more involved (we remark only that in the proof of the recurrence part there 
are some similarities with the proof of Theorem 4.5). 
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5 Frog model on finite graphs and other modifications 

5.1 Finite graphs 

Here we describe some (current) results of the work [4], 

In this section we suppose that Cj is a finite graph, and that the process starts with one-particle-per- 
site initial configuration. In addition, we assume that active particles live some fixed time t, rather than 
having random geometrically distributed lifetime. Now, let Tg be the smallest value t could be, in order to 
guarantee that with probability at least 1 /2 all vertices of Cj will be visited by active particles. The goal is 
now to estimate Tg for different types of graphs. In the sequel we state some of the (preliminary) results 
of [4], 

Theorem 5.1. Let Cj be the the part of the one-dimensional lattice Z restricted to {0,1,2,...,«}. Then 
Tg = 0(log 2 n). 

Theorem 5.2. Let (f be a complete graph with n vertices. Then Tg = O(logn). 

The following result deals with the situation when it is much more difficult to visit all the vertices: 

Theorem 5.3. Let (f be a complete graph with an extra vertex connected by an extra edge to each one of 
its vertices. We have in this situation that Tg = 0(nlog 2 n). 

Next, we present a result (which we hope to improve) for a two-dimensional box: 

Theorem 5.4. Let Cj be the graph with the vertex set {x = {x\ .xf) G Z 2 : 1 < x\ ,X 2 < n} and the edge set 
{[x.y] : \\x — y|| = 1 andx,y £ Z 2 }. Then, Tg is such that 

Clogn < Tg < log 2 n. 

In order to approach the criticality in this model with M = 1, we can perform the above construction 
with FM(Z 2 ,p, 1) instead of FM(Z 2 ,1,1) (i.e., each active particle is killed with probability (1 — p) on 
each step, and, in addition, it is killed when it comes to a site where there is no sleeping particles, see 
Figure 2). 

Open Problem 5.2. Prove that this process survives with positive probability for large enough p. 
Simulations indicate that the critical value of p is around 0.95. 

Now, we describe the second modification. The following model was proposed by Herve Guiol (2000, 
private communication). Consider the frog model in Z 2 , with one-particle-per-site initial configuration, 
but with one important modification of the dynamics of the process: instead of dying, active particles 
become sleeping again with probability (1 — p), on each step (perhaps it is even better to consider this 
model with continuous time, to avoid defining voluntarily what happens in situations when e.g. a particle 
goes asleep in some site, while at the same moment another active particle comes to that site). Then, of 
course, if the original frog model survives for some p, then so does this modified model, but the situation 
with extinction is not at all clear: 

Open Problem 5.3. Prove that this modified frog model becomes extinct for small enough p. 

As we have seen above, the frog model with one-particle-per-site initial configuration always dies out 
in dimension 1 for any p < 1, but here the situation is again not clear: 

Open Problem 5.4. Does the modified frog model survive in dimension 1 with p close enough to 1? 



Frogs and some other interacting random walks models 


285 



Fig. 1 : The modified frog model with M = 1 (and p = 1), after 150 steps. Note that the holes inside the shape will 
never be visited. 
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5.2 Some other modifications 

In this section we describe two natural modifications of FM (Z d ,p, 1), which give rise to some challenging 
open problems. 

Before formulating the first modification, let us recall the shape result from [1] for FM(Z d , 1,1) (cf. 
Section 2). As we noted there, it is generally not possible to identify the asymptotic shape, so simulation 
seems to be the only way to get an idea about it. Performing the simulation for this model is, however, not 
an easy task, because, due to the linear growth speed, at time t there are roughly t d moving particles in 
the system. To simplify the computations, one may note that any particle only encounters a finite number 
of sleeping particles during its lifetime (because an individual random walk has diffusive behavior, while 
the set of visited sites for the frog model grows linearly). So, if during the simulation an active particle 
does not awake anybody for a long time, it is reasonable to suppose that it will not reach the border of the 
growing shape anymore, and so we may remove it from the system without affecting the shape. Thus, let 
us introduce an integer parameter M > 1, and kill any active particle which did not wake up anybody forM 
consecutive times. This idea was recently implemented (in dimension 2) by Fabio Machado and Lucas 
Meyer, and they found that indeed it permitted to speed up the computations. Then, just for curiosity, 
they put M = 1, and the result was quite strange and interesting (see Figure 1): the process survived with 
high probability (at least 0.9), but, apparently, there was no linear growth, and the set of visited sites was 
fractal-like. Now, it is easy to prove that forM large enough the process survives with positive probability, 
but the following problem remains open: 

Open Problem 5.1. Prove that the process with M = 1 survives with positive probability. 



Fig. 2: The modified frog model with M = 1 and p = 0.96, after 300 steps. There are active particles, but only in the 
NW end of the shape. 










Frogs and some other interacting random walks models 


287 


6 Other models with interacting random walks 

In this section we consider the models with Cf = Z‘ I , with infinite lifetime, and where all the particles are 
active from the beginning. 

6.1 Shape theorem 

Let us start from some constant-density initial configuration of particles (it may be e.g. one-particle-per- 
site, or a random configuration governed by Poisson product measure), and let those particles perform 
discrete-time SRWs independently. Moreover, suppose that the particles that are initially in 0 are infected. 
Infected particles transmit the disease to all the particles they meet, and there is no recovering. Similarly 
to Section 2, denote by the set of sites visited by infected particles up to the moment n. The following 
problem seems to be challenging: 

Open Problem 6.1. Prove the analogue of Theorem 2.1 for this model. 

Let us say a few words about the difficulties that arise when one tries to solve the above problem. A 
natural approach would be to consider the time T (0, nx) when the infection comes to the site nx, and then, 
by using the subadditive ergodic theorem, to prove that T {{) .nx) j n converges to some limit as n —> <=°. In 
[1, 3], that was done in the following way: first, a set of auxiliary equally distributed random variables 
T ( kx, (k + 1 )x), k = 0,1,... was constructed in such a way that that sequence is mixing and the subadditive 
inequality 

T(0,nx) < T(0,x) + T(x,2x) 4-1 -T((n — 1 )x,nx) 

is verified for all the realizations of the process. Then, for the random variable T(O.x) it was verified 
that its expectation is finite (that was the difficult part), which permitted us to use the subadditive ergodic 
theorem. Now, for this infection spreading model, it seems to be possible to prove, by using the method 
of [1, 3], that the mean time until the infection comes to x is finite (with a little more work that would 
imply that q„ grows linearly). What is difficult, is to construct a sequence of random variables as above, 
to which some subadditive ergodic theorem is applicable. To clarify a bit more this point, note that, 
even if we start from Poisson product measure (which is invariant with respect to this independent SRWs 
dynamics), the measure at the (random) moment when the infections gets to x for the first time, is not the 
Poisson product measure. 

6.2 A competition model 

Here, exceptionally, we consider simple random walks with continuous time, and the dynamics is de¬ 
scribed as follows: there are red and blue particles; when a red particle comes to a site that contains blue 
particles, all those particles become red, and when a blue particle comes to a site that contains red par¬ 
ticles, all those become blue. Initially, there is one red particle in the origin of Z d , d > 3, and one blue 
particle in each site of an infinite set S cZ‘i The following result was proved in [6]: 

Theorem 6.1. Let x be the random moment when the last red particle becomes blue. We have that X < °° 
a-s. WLxsS ||x|| _(d_2) =°°. 
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